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We investigate the distributions of the link overlap, P{Q), in 3-dimensional Ising spin glasses. 
We use clustering methodology to identify a set of pairs of states from different Gibbs states, and 
calculate its contribution to P(Q). We find that the distribution over this set does not become 
trivial as the system size increases. 



I. INTRODUCTION 

While equilibrium properties of infinite range spin 
glasses^ are completely understood within the framework 
of replica symmetry breaking (RSB)^i2ii^, spin glasses 
with short range interactions are the subject of consid- 
erable current debate and controversy. Open questions 
address the nature of the low temperature phases^i^i^i^ 
and their theoretical description. Resolution of these is- 
sues by experiments or simulations is hindered by the 
extremely long relaxation time required for equilibration. 

The the Edwards- Anderson model is the most widely 
studied type of short-range Ising spin glass 



H{s) 



(ij) 



(1) 



where the sum is over nearest neighbor sites, of a 

simple (hyper) cubic lattice with periodic boundary con- 
ditions. Si — ±1, and the couplings, Jy, are independent 
random variables taken from a normal distribution with 
zero average and standard deviation J = 1. For the in- 
finite range Sherington Kirkpatrick (SK) model the sum 
in Eq. [1] is over all pairs of sites and J = 1/ \/N. 

The system (in 3 or more dimensions) has a finite 
critical temperature. Recent numerical analysis of 3- 
dimensional Ising spin glasses (3DISG) yielded^ Tc ~ 
0.951(9), whereas for the SK model = 1. The high 
temperature phase of the model is a disordered param- 
agnet. As the temperature decreases below Tc the system 
undergoes a transition into a frozen spin-glass phase. 

In the spin glass phase, the microstates are divided 
into Gibbs states; each constitutes an ergodic subset of 
phase space, i.e. a maximal subspace that the system 
can span (or visit) as time tends to infinity. In a finite 
system phase space consists of one such state; however, 
we identify the infinite volume Gibbs states with subsets 
of the phase space surrounded by free energy barriers, 
whose height diverges as iV ^ oo. Here the term "Gibbs 
states" refers to such subsets. 

Extensive recent numerical investigations, done at fi- 
nite riSiiiii^, as well as ground state analysisi^^ii^ii^, 
suggest evidence for a multiplicity of Gibbs states in the 
low temperature phase of the 3DISG. The most widely 



measured properties are 
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the site overlap between any two microstates s and t, the 
global distribution of q, 

pM = ^ e-«(«) ^-"''^5{q{s, t)-q), (3) 



and, in particular, p{q) = pj{q), the distribution aver- 
aged over the random variables {Jij}- In ^ Z is the 
partition function of the system. 

There is general agreement that in 3DISG the averaged 
distribution p((7) is not trivial (that is, does not converge 
to a (5-function in the infinite volume limit) and is not 
self averaging, in agreement with RSB theory. 

As to the droplet theorj*^:^, although its common in- 
terpretation involves a trivial p{q)^, it explicitly predicts 
only the triviality of the local p(q) - the overlap distri- 
bution over a finite (large) window in an infinite system. 
That is, locally there is only one Gibbs state for the sys- 
tem (up to a spin- flip symmetry), so that when a finite 
window of an infinite system is viewed, the system will 
be almost always in this Gibbs state. 

In order to test this prediction numerically, one should 
observe the site overlap in a finite constant window as 
the system size increases^*^. An alternative is to measure 
the link overlap 
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summing over all the iVb bonds in the system, e.g. over 
all nearest neighbors pairs in the case of 3DISG. The 
distribution of the link overlap, Pj((3) = {S{Q{s,t) — Q)), 
is defined similarly to p.j{q) in Eq. [31 The average over 
realizations is P{Q) — Pj{Q). 

According to the droplet theory P{Q) is trivial and 
consequently Pj{Q) — P{Q)- This was predicted earlier 
from scaling analysis of numerical results22i2i. Newman 
and Stein^i^^ showed that triviality of P{Q) for given 
boundary conditions may be deduced from general con- 
siderations even if one relaxes some of the scaling as- 
sumptions of the droplet theory. However, according to 
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FIG. 1: The variance of the hnk overlap distribution, P{Q), 
for 3DISG systems. The error bars are larger then the values 
of the variance. 



RSB P{Q) is not trivial. Unlike the case of p{q), the 
triviality of P{Q) is still an open questiort^ i^'^i^^i^^ . 



II. NUMERICAL SIMULATIONS 

In the present work we analyze 3DISG systems of sizes 
L = 8, 10, 12. We generated 65 realizations of the dis- 
order, {Jy}, for each system size. For each realization 
we produced a weighted sample of microstates. We used 
simulated tempering Monte Carlo^^. The number of tem- 
peratures, Nt, and their values were determined so that 
in each inter-temperature step the transition probability 
is larger than 0.25. 

We used Nt = 16 for i = 8, and Nt > 27 for L = 
10, 12. The lowest and highest temperatures were Ti = 
0.24J/A:b and Tnj, = 2J/kB {Tn^ = 0.92J/kB for L = 
8). For every realization, we took TOOiVr samples (i.e. 
microstate-temperature pairs). We ran 20Nt^ sweeps to 
thermalize the system, and the same number of sweeps 
between samples. 

In most realizations the temperatures of consecutive 
samples were decorrelated, so that the system often went 
above Tc, where it lost spin correlations. For some real- 
izations the decorrelation times were as large as 10 sam- 
pling cycles. 



III. LINK OVERLAP DISTRIBUTION 

One approach used to study the link overlap is to mea- 
sure the fractal dimension of the surface of low energy 
excitationsiiiSi. Such studies are limited by the small 
system sizes available for numerical investigation. 

Another approach is to calculate directly the moments 
of the link overlap distribution, P(Q)^^-^*. In Fig. [l]we 

present (Q^) — (Q) , the second moment of P{Q)- As 
the system size increases, the second moment decreases 
(apparently towards zero) at all temperatures, indicating 
that P{Q) becomes trivial a.s L ^ oo. 

However, characterization of the behavior of a distribu- 
tion on the basis of a single parameter, such as a fractal 
dimension or a particular moment, may be misleading. 



This was demonstrated in Ref. |29| for the site overlap 
distribution, and the same holds for the link overlap dis- 
tribution. 

The distribution Pj{Q), measured for a particular 
realization of the random {Jij}, may consist of many 
modes, which correspond to different peaks in the graph 
of Pj{Q)- Each such mode is related to the distribution 
of overlaps between two given Gibbs states. We denote 
by Gi the Gibbs state with the highest statistical weight. 

As temperature decreases the variation of the statis- 
tical weight of the Gibbs states increases. In particular, 
the weight of Gi approaches 1 as T goes to zero. Thus, 
at low temperatures the most dominant mode of Pj{Q) 
is 

Pi(Q) =Prob(Q(s,t) =Q|s,t e Gi) (5) 

- the distribution of overlaps within Gi - which is triv- 
ial in the thermodynamic limit. At small system sizes, 
as the ones considered in this work, the peak of Pi{Q) 
is still wide, and it narrows as L increases. Since even 
at finite (albeit low) temperatures Pi{Q) dominates the 
distribution, one would expect a clear decrease of the 
variance of P{Q) with size, as observed in Fig. [1] even if 
P{Q) is non-trivial in the thermodynamic limit. Hence 
decrease of the second moment of the distribution with 
size may capture nothing but the narrowing of the dom- 
inant single peak part of a non-trivial P{Q)', in order to 
test non-triviality, it is much more effective to measure 
the size dependence of a quantity which is independent 
of the dominant (but uninformative) Pi{Q). 

We present P{Q) in Fig. Efa) for two temperatures, 
one above and one below Tc, for three system sizes. 
Above Tc, P{Q) consists of a single peak which narrows 
as the L increases;, approaching 5{Q — {Q)). Below Tc 
the distribution consists of a dominant peak, whose main 
contribution comes from Pi(Q), and a tail at lower val- 
ues of Q. As L increases, the peak narrows and since the 
weight of the peak is much larger then that of the tail, 
a single measured parameter (say moment) may reflect 
mainly the convergence of the peak to a (5 function, which 
is expected irrespective of whether P{Q) is trivial or not. 
To address this question one has to assess whether the 
weight of the tail of P{Q), say, in the range < Q < 0.8, 
decreases with increasing system size. 

If P{Q) is to have a non-trivial component at the ther- 
modynamic limit, it must originate from pairs of mi- 
crostates that belong to different Gibbs states. Ideally, 
we would like to calculate the link overlap distribution 
over all such pairs. In finite systems this is not possible 
since we cannot identify the Gibbs states unambiguously. 
However, the coarsest partitions of phase space can be 
observed clearly, even for very small systems^^'^^. This 
partition consists of two sets of microstates, Ci and C2, 
which we identify by the average linkage agglomerative 
clustering algorithm^. 

The statistical weights of the sets Ci and C2 remain 
finite as the system size increases, and the average site 
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FIG. 2: (a) The link overlap distribution, P(Q) measured 
for 3DISG systems at T = 1.35J/fcB > Tc The distribution 
converges to a 5-peak. The grey lines above and below each 
curve indicate the error range. 



overlap between them remains different from the over- 
lap between two microstates that belong to the same 
se^. Therefore the contribution of the set of pairs of 
microstates {{s,t) \s £ Ci,t G C2} to P{Q) remains fi- 
nite as the system size increases. 

Following the method described in Ref. 29 for p{q), we 
define 

Po(Q) - Prob( Q(s, t) = g I s e Ci, t e C2 ) . (6) 

In order to move from a state in Ci to a state in C2, 
some groups of spins (domains) have to be flipped. We 
refer to the surfaces that separate these flipped domains 
from spins that did not flip as domain walls. If these 
domain walls have a vanishing density, the fraction of 
links affected by switching one micro-state from Ci to C2 
goes to zero as the system size increases, and hence at 
the thermodynamic limit we will have Po{Q) — P{Q)- 

In Fig. [Sfa) we observe that for 3DISG at T = 
0.28J//cb, P{Q) seems to converge to a peak at Q ~ 0.9. 
If Po{Q) does not converge to the same distribution, we 
expect it to have a non-vanishing weight at Q < 0.9. 
We calculate /°-® Po{Q)dQ. For L = 8, 10, 12 the values 
of the integral are 0.81(5), 0.83(4), 0.80(4) respectively, 
suggesting that the weight of PoiQ) for < Q < 0.8 does 
not decrease. Thus, Po{Q) does not converge to the peak 
at Q ~ 0.9 which dominates P{Q)- 

These results imply that the domain walls between Ci 
and C2 occupy a finite fraction of the volume of the sys- 
tem. Since Ci and C2 remain finite and distinguishable as 
the system size increasesi^, Po{Q) (shown in Fig. El^b)) 
constitutes a finite contribution to P{Q). Consequently, 
P{Q) is not trivial at the thermodynamic limit. 




FIG. 3: (a) The link overlap distribution measured for 3DISG 
systems at T = 0.28J/feB < 7c, in linear plot and linear-log 
plot (inset), (b) The distribution Po{Q) for T = 0.28J/fcB, 
which gives a finite contribution to P(Q)- All distributions 
are normalized to unity. The grey lines above and below each 
curve indicate the error range. 



IV. DISCUSSION 



According to RSB theory, P{Q) is not trivial. More 
over, all overlap measures are equivalents^!^, so Q is de- 
termined by the value of the site overlap q. This state- 
ment is trivial for the SK model, where the sum in Eq. [4] 
is over all pairs of sites, soQ{s,t) = {q^{s, t) — N~^)/{\ — 
N^^). In 3DISG the site overlap and bond overlap have 
different roles: given two states of the system, the site 
overlap is given by the volume of spin-domains flipped 
between the states, and the link overlap by the surface 
area of the domain walls. 

Recent numerical investigation of 3DISG systems2^ in- 
dicated that the two overlaps are equivalent. This re- 
sult along with previous numerical evidence for the non- 
triviality of the site overlap, p{q), yield a non-trivial P{Q) 
and support the results presented here. 
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